THE KISSING NUMBER IN FOUR 
DIMENSIONS 

Oleg R. Musin * 

Abstract 

The kissing number problem asks for the maximal number k{n) of 
equal size nonoverlapping spheres in n- dimensional space that can touch 
another sphere of the same size. This problem in dimension three was the 
subject of a famous discussion between Isaac Newton and David Gregory 
in 1694. In three dimensions the problem was finally solved only in 1953 
by Schiitte and van der Waerden. 

In this paper we present a solution of a long-standing problem about 
the kissing number in four dimensions. Namely, the equality fc(4) — 24 is 
proved. The proof is based on a modification of Delsarte's method. 

1 Introduction 

The kissing number kin) is the highest number of equal nonoverlapping spheres 
in R" that can touch another sphere of the same size. In three dimensions the 
kissing number problem is asking how many white billiard balls can kiss (touch) 
a black ball. 

The most symmetrical configuration, 12 billiard balls around another, is if 
the 12 balls are placed at positions corresponding to the vertices of a regular 
icosahedron concentric with the central ball. However, these 12 outer balls do 
not kiss each other and may all moved freely. So perhaps if you moved all of 
them to one side a 13th ball would possibly fit in? 

This problem was the subject of a famous discussion between Isaac Newton 
and David Gregory in 1694. It is commonly said that Newton believed the 
answer was 12 balls, while Gregory thought that 13 might be possible. However, 
Casselman (8j found some puzzling features in this story. 

The Newton- Gregory problem is often called the thirteen spheres problem. 
Hoppe JS] thought he had solved the problem in 1874. However, there was 
a mistake - an analysis of this mistake was published by Hales jlTj in 1994. 
Finally, this problem was solved by Schiitte and van der Waerden in 1953 |31| . 
A subsequent two-page sketch of a proof was given by Leech W3 in 1956. The 
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thirteen spheres problem continues to be of interest, and several new proofs have 
been published in the last few years pOlOl BllTldB ] . 

Note that fc(4) > 24. Indeed, the unit sphere in R'* centered at (0,0,0,0) 
has 24 unit spheres around it, centered at the points (±-\/2, ±\/2, 0, 0), with any 
choice of signs and any ordering of the coordinates. The convex hull of these 24 
points yields a famous 4-dimensional regular polytope - the "24-ceH" . Its facets 
are 24 regular octahedra. 

Coxeter proposed upper bounds on k{n) in 1963 jlQi; for n = 4,5,6, 7, and 
8 these bounds were 26, 48, 85, 146, and 244, respectively. Coxeter's bounds 
are based on the conjecture that equal size spherical caps on a sphere can be 
packed no denser than packing where the Delaunay triangulation with vertices 
at the centers of caps consists of regular simplices. This conjecture has been 
proved by Boroczky in 1978 [S]. 

The main progress in the kissing number problem in high dimensions was 
made in the end of 1970s. In 1978: Kabatiansky and Levenshtein have found 
an asymptotic upper bound 2° ''°^"(^+°(i)) for fc(n) |23. (Currently known the 
lower bound is 20-2075n(H-o(i)) |32j ) 1979. Levenshtein [331, and indepen- 
dently Odlyzko and Sloane [221 (= O Chap. 13]), using Delsarte's method, have 
proved that fc(8) = 240, and k{2A) = 196560. This proof is surprisingly short, 
clean, and technically easier than all proofs in three dimensions. 

However, n = 8, 24 are the only dimensions in which this method gives a 
precise result. For other dimensions (for instance, n = 3,4) the upper bounds 
exceed the lower. In j^T] the Delsarte method was applied in dimensions up 
to 24 (see Table 1.5]). For comparison with the values of Coxeter's bounds 
on fc(n) for n = 4,5,6,7, and 8 this method gives 25, 46, 82, 140, and 240, 
respectively. (For n ~ 3 Coxeter's and Delsarte's methods only gave fc(3) < 13 

[101123.) 

Improvements in the upper bounds on kissing numbers (for n < 24) were 
rather weak during next years (see |5| Preface to Third Edition] for a brief review 
and references). Arestov and Babenko proved that the bound fc(4) < 25 
cannot be improved using Delsarte's method. Hsiang ^\ claims a proof of 
fc(4) = 24. His work has not received yet a positive peer review. 

If M unit spheres kiss the unit sphere in R", then the set of kissing points 
is an arrangement on the central sphere such that the (Euclidean) distance 
between any two points is at least 1. So the kissing number problem can be 
stated in other way: How many points can be placed on the surface of S"^^ so 
that the angular separation between any two points is at least 7r/3? 

This leads to an important generalization: a finite subset X of S"~^ is called 
a spherical ip-code if for every pair {x,y) of X the inner product x ■ y < cosip, 
i.e. the minimal angular separation is at least ip. Spherical codes have many 
applications. The main application outside mathematics is in the design of 
signals for data transmission and storage. There are interesting applications to 
the numerical evaluation of n-dimensional integrals 9, Chap. 3]. 

Delsarte's method (also known in coding theory as Delsarte's linear program- 
ming method or Delsarte's scheme) is widely used for finding bounds for codes. 
This method is described in E] (see also |2H1 for a beautiful exposition) . 
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In this paper we present an extension of the Delsarte method that allowed 
to prove the bound fc(4) < 25, i.e. fc(4) — 24. This extension yields also a proof 
for fc(3) < 13 Uni- 

The first version of these proofs used numerical solutions of some nonconvex 
constrained optimization problems [23] (see also (2H|)- Now, using geometric 
approach, we reduced it to relatively simple computations. 

The paper is organized as follows: Section 2 shows that the main thorem: 
fc(4) = 24 easily follows from two lemmas: Lemma A and Lemma B. Section 3 
reviews the Delsarte method and gives a proof of Lemma A. Section 4 extends 
Delsarte's bounds and reduces the upper bound problem for i/)-codes to some 
optimization problem. Section 5 reduces the dimension of the corresponding 
optimization problem. Section 6 develops a numerical method for a solution of 
this optimization problem and gives a proof of Lemma B. 

Acknowledgment. 1 wish to thank Eiichi Bannai, Dmitry Leshchiner, Sergei 
Ovchinnikov, Makoto Tagami, Giinter Ziegier, and especially anonymous refer- 
ees of this paper for helpful discussions and useful comments. 

1 am very grateful to Ivan Dynnikov who pointed out a gap in arguments on 
earlier draft of [25] . 

2 The Main Theorem 

Let us introduce the following polynomial of degree nine:^ 

1344 „ 2688 , 1764 , 2048 , 1229 . 516 . 217 2 

fJt) := f H H f t 

■' ^ ' 25 25 25 125 125 125 500 125 

Lemma A. Let X = {xi, . . . , xm} be points in the unit sphere . Then 

M M 
i=l = 1 

We give a proof of Lemma A in the next section. 

Lemma B. Suppose X = {xi, . . . ,xm} is a subset of such that the angular 
separation between any two distinct points at least 7r/3. Then 

M M 

5(X)-^^/4(a;, •x,)<25M. 

i=l 3 = 1 

A proof of Lemma B is given in the end of Section 6. 
Main Theorem. fc(4) = 24. 

^Tho polynomial was found by the linear programming method (see details in the 
Appendix). This method for n = 4, ^ = 1/2, d = 9, N = 2000, to = 0.6058 gives E x 24.7895. 
For fi coefficients were changed to "better looking" ones with E 24.8644. 
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Proof. Let X be a spherical 7r/3-code in with M = fc(4) points. Then X 
satisfies the assumptions in Lemmas A, B. Therefore, < S{X) < 25M. 
From this M < 25 follows, i.e. M < 24. From the other side we have A:(4) > 24, 
showing that M = fc(4) = 24. □ 



3 Delsarte's method 

From here on we will speak of a; G S"""'^ alternatively of points in S"""'^ or of 
vectors in R". 

Let X = {xi,X2, ■ ■ ■ , xm} be any finite subset of the unit sphere S"^^ C R", 
S"~-^ — {x : X £ R", X ■ X — = 1}. By (f)ij = dist{xi,Xj) we denote the 

spherical (angular) distance between Xi, Xj. Clearly, cos — Xi ■ Xj. 

3-A. Schoenberg's theorem. Let ui, U2, . . . , um be any real numbers. Then 

ll^WiXjlp = cos (t)i J UiU J > 0, 
hi 

or equivalently the Gram matrix (^cos (jjij^ is positive semidefinite. 

Schoenberg [221 extended this property to Gegenbauer polynomials G^"^ . He 
proved that the matrix ^G^""* (cos 0^ ,,)^ is positive semidefinite for any finite 

X C S"-i. 

Schoenberg proved also that the converse holds: if f{t) is a real polynomial 
and for any finite X C S"^^ the matrix (/(cos j)) is positive semidefinite, 

then f{t) is a linear combination of G''l^\t) with nonnegative coefficients. 

3-B. The Gegenbauer polynomials. Let us recall definitions of Gegenbauer 
polynomials. Let polynomials cj^^^t) are defined by the expansion 

oc 

(l_2rt + r2)(2-«)/2^^r'=ci")(i). 

fc=0 

Then the polynomials G^^^t) :— G^"'''(t)/G|,"''(l) are called Gegenbauer or ul- 
traspherical polynomials. (So the normalization of g[,"^ is determined by the 
condition g[,"^(1) = 1.) Also the Gegenbauer polynomials G^"^ can be defined 
by the recurrence formula: 

/-.(")_ 1 M"-) - ^- ^(n) _ + - 4)t G^"\ - (fc - 1) G^"^2 

Go -1, Gi G, - fc + „_3 

They are orthogonal on the interval [—1,1] with respect to the weight func- 
tion p{t) = (1 - t2)(n-3)/2 (ggg details in |7||ni[TniEni)- In the case n = 3, Gf' 
are Legendre polynomials Ph , and G^'''' are Chebyshev polynomials of the second 
kind (but with a different normalization than usual, Vk{X) — 1), 
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For instance, Uq = 1, Ui = t, U2 ^ (U"^ - l)/3, U3 = 2t^ - t, 
Ui = (let'' - 12^2 + l)/5, . . . , C/g = (256f9 - 512f + SSGt'^ - 80t^ + 5i)/5. 

3-C. Delsarte's inequality. If a symmetric matrix A is positive semidefinite, 
then the sum of all its entries is nonnegative. Schoenberg's theorem implies that 
the matrix is positive semidefinite, where tij := cos, (f)ij, Then 

M M 

Y.llGt\ti,)>^ (3-1) 

i=l J = l 

Definition 1. We denote by t/ie set of continuous functions f : [— 1, 1] — + R 
representable as series 

00 

f{t)=J2ckGi"\t) 

fe=0 

whose coefficients satisfy the following conditions: 



Co > 0, Ck> for k = 1,2,..., /(I) = J^Cfc < oo. 

fe=o 

Suppose / e G+ and let 

M M 

Using (3.1), we get 

00 / M M \ M M 

S{X) = ^cJj2J2Gt'\hj)\>E12^Go"\tij) = coM'. 

k=0 \i=l j=l J i=l j=l 

Then 

S{X) > coM^. (3.2) 

3-D. Proof of Lemma A. 

Proof. The expansion of fi in terms of Uk = G^^^ is 

153,, 871,, 128,, 21 , 
f. = U, + 2U,^-U, + -U, + -U,^-U, 

We see that f^ e G4 with cq = 1. So Lemma A follows from (3.2). □ 

3-E. Delsarte's bound. Let X = {xi, . . . ,xm} C S"""' be a spherical V'-codc, 
i.e. for all i 7^ j, tij = cos^ij = Xi ■ Xj < z := cos^', i.e. tij G [—^,z] (but 
ti,i = l). 
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Suppose / e G+ and f{t) < for all t e [-1,^], then /(^ij) < for all 
i =/= j. That implies 

Sf{X) = Mf{l) + 2/(ti,2) + . . . + 2f{tM-i,M) < Mf{l). 

If we combine this with (3.2), then we get M < /(l)/co. 

Let A(n,'0) be the maximal size of a tp-code in S"~^. Then we have: 

ffl) 

A(n»<^ (3.3) 
Co 

The inequality (3.3) play a crucial role in the Delsarte method (see details 
in El mini ESI m Ell ESI EH )■ W ^ = 1/2 and CO = 1, then (3.3) implies 

k{n) = yl(n,7r/3) < /(I). 

Levenshtein j^Sl, and independently Odlyzko and Sloane ITTI for n = 8, 24 have 
found suitable polynomials f{t): f{t) < for all t e [—1, 1/2], / e G+, cq = 1 
with 

/(I) = 240 for n = 8; and /(I) = 196560 for n = 24. 

Then 

fc(8) < 240, A:(24) < 196560. 

For n = 8, 24 the minimal vectors in sphere packings Eg, and Leech lattice give 
these kissing numbers. Thus k{8) = 240, and A:(24) = 196560. 

When n = 4, a polynomial / of degree 9 with /(I) « 25.5585 was found in 
123. This implies 24 < fc(4) < 25. 

4 An extension of Delsarte's method. 

4-A. An extension of Delsarte's bound. 

Let f{t) be any real function on the interval [—1,1]. Let for a given ip 
z := cos V". Consider on sphere S"~^ points yo, yi, . . . , ym such that 

Vi -Vj < z for all i ^ j, f{yo ■ yi) > for 1 < i < m. (4.1) 

Definition 2. For fixed yo £ S"^^, m > 0, z, and f{t) let us define the family 
Qm{yo) — Qm{yo,n-, f) of finite sets of points from S""^ by the formula 

n (1, ■) — / 'fyoj'' TO = 0, 

<^M) - I |y ^ ^ j^^j ^ s«-i : {yo} U Y satisfies (4.1)}, to > 1. 

Denote fi = ^(n, z, f) := max{m : Qm{yo) 7^ 0}- 

For Q < m < fi we define the function H — Hf on the family Qmiyo)-' 

H{yo):=f{l) for to = 0, 
H{yo;Y) = H{yo;yi, ...,?/„) := /(I) + f{yQ ■ y^) + . . . + f{yo ■ y™) for to > 1. 
Let 

h„i^hm{n,z,f):= sup {H{yo;Y)}, hmax ■= max{ho,hi, . . . ,hfj,}. 
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Theorem 1. Suppose / e G+. Then 

I. ^ hmax{n,C0S1p,f) 1 

< = — max|/io, hi,. .. , ha\- 

Co Co 

Proof. Let X = {a;i, . . . ^Xm} C S"~^ be a spherical V-code. Since / € G+, 
(3.2) yields: S{X) > cqM"^. 

Denote J{i) := {j : f{xi ■ xj) > 0, j ^ i}, X{i) := {xj : j G J{i)}. Then 

M 

Si{X) := J2 fi^i ■ ^j) < /(I) + Yl • ^^•) = -^(^^5 ^(0) < ^ma., 

J=i jeJii) 

so then 

M 

= J2 Si{X) < Mhma.. (4.2) 

i=l 

We have CqM^ < S{X) < Mhmax, i-C cqM < /iman as required. □ 

Note that ho = /(I). If f{t) < for ah t G [-1,^], then iJL{n,z,f ) = 0, 
i.e. hmax = ho = /(I). Therefore, this theorem yields the Delsarte bound 
M < /(l)/co. 

4-B. The class of functions <P(to,z). 

The problem of evaluating of hmax in general case looks even more compli- 
cated than the upper bound problem for spherical V'-codes. It is not clear how 
to find /i, what is an optimal arrangement for Y7 Here we consider this problem 
only for a very restrictive class of functions <l>{to,z). For the bound given by 
Theorem 1 we need / G G+. However, for evaluations of hm we don't need this 
assumption. So we are not assume that / G G+. 

Definition 3. Let real numbers to,z satisfy 1 > to > z > Q. We denote by 
^{to, z) the set of functions f : [—1, 1] — » R such that 

f{t)<0 for te[-to,z]. 

Let / G ^{to, z), and let Y G flrnivo^n; /). Denote 

eo := — yoj ^0 arccosto, ^'i dist(eo, y^) for i=l,...,m. 

(In other words, cq is the antipodal point to yo-) 

It is easy to see that /(yo • 2/i) > only if 6i < 9o- Therefore, 

Y is a spherical tp-code in the open spherical cap Cap(eo,^o) of center eo and 
radius 9o with Tr/2 > ip > 6o. 

This assumption is quit restrictive and in particular derives the convexity prop- 
erty for Y. We are using this property in the next section. 

4-C. Convexity property. A subset of S""-'^ is called spherically convex if it 
contains, with every two nonantipodal points, the small arc of the great circle 



7 



containing them. The closure of a convex set is convex and is the intersection 
of closed hemispheres (see details in |12l). 

Let Y — {j/i, . . . , Urn} C Cap(eo, 9q), Oq < tt/2. Then the convex hull of Y 
is well defined, and is the intersection of all convex sets containing Y. Denote 
the convex hull of Y by = Am{Y). 

Recall a definition of a vertex of a convex set: A point y £ W is called the 
vertex (extremal point) of a spherically convex closed set W, if the set W \ {y} 
is spherically convex or, equivalently, there are no points x, z from W for which 
y is an interior point of the minor arc xz of large radius connecting x, z. 

Theorem 2. Let Y — {yi, . . . , C S"^^ be a spherical %jj-code. Suppose 
Y C Cap(eo,0o), and < 6q < ijj < 7r/2. Then any y^. is a vertex 0/ A,„. 

Proof. The cases m = 1,2 are evident. For the case to = 3 the theorem can be 
easily proved by contradiction. Indeed, suppose that some point, for instance, 
2/2 is not a vertex of A3. Then, firstly, the set A3 is the arc yiys, and, secondly, 
the point y2 lies on the arc yiys- From this it follows that dist(?/i, j/3) > 2ip, 
since y is a ■0-code. From the other hand, according to the triangle inequality, 
we have 

2^ < dist(yi,y3) < dist(eo,yi) +dist(eo,y3) < 26*0. 

We obtained the contradiction. It remains to prove the theorem for to > 4. 

In this paper we need only one fact from spherical trigonometry, namely the 
law of cosines (or the cosine theorem): 

cos (j) = cos 9i cos 62 + sin 61 sin 6*2 cos (p, 

where for a spherical triangle ABC the angular lengths of its sides are 
dist(A, B) = 9i, dist(A, 0)^02, dist(B, C) = (/), and ZBAC = ip. 

By the assumptions: 

Ok = dist(j/fe, eo) <9a <-tJj for 1 < fc < to; (/i^j := dist(yfc, y^) > -0, ^ 7^ j- 

Let us prove that there is no point belonging both to the interior of Am and 
relative interior of some facet of dimension c?, I < d < dim Am ■ Assume the 
converse. Then consider the great (n — 2)-sphere flk such that yk G ^k, and 
rjfc is orthogonal to the arc eoyk- (Note that 9k > 0. Conversely, yk = cq and 
(j)k,j = 9j <9o < ip.) 

The great sphere ftk divides S"^^ into two closed hemispheres: Hi and H2. 
Suppose Co lies in the interior of Hi , then at least one yj belongs H2 ■ Consider 
the triangle eaykyj and denote by "fkj the angle ZeoykUj in this triangle. The 
law of cosines yields 

cos 9j = cos 9k cos (jjkj + sin 9k sin (pkj cos jk,j 

Since yj £ H2, we have jkj > 90°, and cosjkj < (Fig. 1). From the 
conditions of Theorem 2 there follow the inequalities 

sm9k > 0, sin(j)kj > 0, cos9k > 0, cos9j > 0. 
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Hence, using the cosine theorem we obtain 

cos 9j = cos 9k cos (pkj + sin 9k sin (pkj cos 'jkj , 

< cos 9j < cos 9k cos (()k,j ■ 
Prom these inequaUties and < cos^fc < 1 there follow that, firstly, 

< cos(j)kj ^i.e. ip < 4>k,j < t/^)) 
and, secondly, the inequalities 

COS^j < cos (l)k,j < cosip. 
Therefore, 6j > tp- This contradiction completes the proof of Theorem 2. □ 

eo 




Fig. 1 




4-D. Bounds on /i. 

Theorem 3. Let Y = {yi, . . . ,ym} C S""^ be a spherical ^-code. Suppose 
Y c Cap(eo, ^o), and < ip/2 <9o<ip< 7r/2. Then 



m < A [ n — 1, arccos 



cos -0 — COS^ 00 
sin^ ^0 



Proof. It is easy to see that the assumption < tp/2 < 9o < tp < Tr/2 
guarantees, firstly, that the right side of the inequality in Theorem 3 is well 
defined, secondly, that there is Y with to > 2. 

If TO > 2, then yi ^ eo- Conversely, ip < dist{y i,yj) — dist(eo, yj) — 9j < 9o, 
a contradiction. Therefore, the projection 11 from the pole eo which sends 
X G S"~^ along its meridian to the equator of the sphere is defined for all y,. 

Denote jij := dist {U{yi),U{yj)) (see Fig. 2). Then from the law of cosines 
and the inequality cos <pij < z = cos ip, we get 

cos O, , — cos 0, cos 9j Z — cos Oi cos Oj 
COSJiJ = — z z < z z 



sm Ui sm Uj 



sm Ui sm f 
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Let 



^, 2: — cos a cos /3 , dR(a,B) cos/3 — zcosa 
R{a,l3)^ . .„ ^ , then = — ^ 



sin a sin P ' cJa sin a sin /3 

We have 6*0 < ip. Therefore, if < a, f3 < 9o, then cos/3 > z. That yields: 
dR{a, (3) /da > 0, i.e. R{a, /3) is a monotone increasing function in a. We 
obtain Ria,l3) < R{0o,l3) = i?(/3,0o) < RiOo,0o). 
Therefore, 



z — cos t/i cos W, z ~ cos t/p 

cos 7i , < — - < 5 = cos 0. 

sm Oi sm t^j sin Oq 



Thus n(F) is a ^-code on the equator S" ^. That yields m < A{n — 1, (5). □ 

Corollary 1. Suppose f G <P{tQ,z). If 21^ > z + 1, then ij{n,z,f) — 1, 
otherwise ^ 

fi{n,z,f) < A(n— 1 , arccos 1 j , 

Proof. Let cosV' = z, cos 6*0 = to- Then 2tQ > z + 1 if and only if > 26*0. 
Clearly that in this case the size of any 1/1-code in the cap Cap(eo, ^0) is at most 
1. Otherwise, ip < 29q and this corollary follows from Theorem 3. □ 

Corollary 2. Suppose f e (I>{to,z). Then 

Ai(3,z,/) < 5. 

Proof Note that 

T = ^— ^ < ^ ~ = < -. Then 5 = arccosT > tt/3. 

I -tl - I - z^ 1 + z 2 ' 

Thus /i(3, z, /) < A(2, 5) < 2tt/5 < 6. □ 

Corollary 3. Suppose f E (I>{to,z). 

{{) If to > V^, then /i(4, z, /) < 4. 

(m) // z = 1/2, to > 0.6058, t/ieK Ai(4,z,/) < 6. 



Proof. Denote by (pk{M) the largest angular separation that can be attained 
in a spherical code on S*"'"^ containing M points. In three dimensions the best 
codes and the values ip3{M) presently known for Af < 12 and M — 24 (see 
^2 El EDI)- Schiitte and van der Waerden 30 proved that 

(^3(5) = ^3(6) = 90°, cos (^3 (7) = cot 40° cot 80°, (/sg (7) w 77.86954°. 

(i) Since z-t^<0, Corollary 1 yields: ^{4, z, f) < A{3, S), where S > 90°. We 
have S > <P3(5). Thus ^ < 5. 
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(m) Note that for to > 0.6058, 

1 /2 _ f2 
arccos-( ^ > 77.87°. 

So Corollary 1 implies /i(4,l/2,/) < ^(3,77.87°). Since 77.87° > ifsi?), we 
have ^(3, 77.87°) < 7, i.e. ii<6. □ 

4-E. Optimization problem. Let 

to'-=cos6o, z:=cosV', cos(5 := 1, := A(n — 1, i5). 

1 — 

For given n, i/', Oq, f G ^(^o, -z), € S"~-^, and m < fi* , the value hm{n, z, /) is 
the solution of the following optimization problem on S"~^: 

maximize /(I) + /(-eo • yi) + . . . + /(-eo • Vm) 

subject to the constraints 

Vi e S""\ i = 1, . . . , m, dist(eo, Vi) < 6o, dist(?/i, yj) > V, « i- 

The dimension of this problem is (n— l)m < (n— If /x* is small enough, 
then for small n it gets relatively small - dimensional optimization problems 
for computation of values hm- If additionally f{t) is a monotone decreasing 
function on [—1, —to], then in some cases this problem can be reduced to {n—1) 
- dimensional optimization problem of a type that can be treated numerically. 



5 Optimal and irreducible sets 

5-A. The monotonicity assumption and optimal sets. 

Definition 4. We denote by ^*{z) the set of all functions / S U ^{to,z) 

To>Z 

such that f{t) is a monotone decreasing function on the interval [— 1,— tq], 
and /(-I) > > /(-To). 

For any f G $*{z), denote to = to{f) := sup{t e [tq, 1] : f{-t) < 0}. 

Clearly, if / e $*{z), then / e ^{to,z), i.e. /(t) < for t e [-1,-io]- 
Moreover, if /(f) is a continious function on [—1, — z], then f{—to) = 0. 

Consider a spherical V'-code Y — {yi, . . . , ym} C Cap(eo, Oo) C S"^^. Then 
we have the constraint: (pij := dist{yi,yj) > ip for all i ^ j. Denote by T\f,{Y) 
the graph with the set of vertices Y and the set of edges yiyj with = tp. 

Definition 5. Let f € ^*{z), V' = arccos(z), 6^0 = arccos(fo)- We say that 
a spherical ijj-code Y = {yi, . . . ,ym} C Cap(eo,^o) C S"~^ is optimal for f if 
Hf{-eo; Y) = h„^{n, z, /). 

// optimal Y is not unique up to isometry, then we call Y as optimal if the 
graph T^{Y) has the maximal number of edges. 
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Let 6k ■= dist(j/fc, eo). Then H{—eo;Y) can be represented in the form: 



FfiOu ...,9m) := Hf{-eo; Y) = /(I) + /(- cos^i) + ... + /(- cos^„). 



Let us call F{9i, . . . , 0m) = Ff{9i,. . . , 9m) the efficient function. Clearly, if 
/ e ^*{z), then the efficient function is a monotone decreasing function in the 
interval [0,9o] for any variable 9k. 

5-B. Irreducible sets. 

Definition 6. Let < 9o < tp < tt/2. We say that a spherical ip-code Y = 
{yi, . . . , ym} C Cap(eo, 9o) C S"~^ is irreducible (or jammed) if any yk can not 
be shifted towards eo (i.e. this shift decreases 9k) such that Y' , which is obtained 
after this shifting, is also a ip-code. 

As above, in the case when irreducible Y is not defined uniquely up to isom- 
etry by 9i, we say that Y is irreducible if the graph T^(Y) has the maximal 
number of edges. 

Proposition 1. Let f e ^*(z). Suppose Y c Cap(eo,^o) C S"~^ is optimal 

for f . Then Y is irreAucible. 

Proof. The efficient function F{9\, . . . ,9m) increases whenever 9k decreases. 
Prom this follows that yk can not be shifted towards cq. In the converse case, 
Go; Y) ~ F{9i, . . . , 9m) increases whenever yk tends to cq. It contradicts 
the optimality of the initial set Y . □ 

Lemma 1. IfY = {yi, . . . , ym} is irreducible, then 

(i) eo G =convex hull ofY; 

(ii) If m> 1, then degj/j > for all y, e Y, where by degy, denoted the degree 
of the vertex yi in the graph r^(y). 

Proof, (i) Otherwise whole Y can be shifted towards cq. 

(ii) Clearly, if (f>ij > ip for all j ^ i, then yi can be shifted towards eo. □ 

For m = 1 from this follows that eo = yi, i.e. hi = s\x\>{F(9\)} = F(0). 
Thus 



Consider A„ c S"~^ of dimension k, dimA^ = k. Since A^ is a convex 
set, there exists the great A;-dimensional sphere S*^ in S"""'^ containing A^. 

Note that if dim A^ = 1, then m = 2. Indeed, since dimA^ = 1, it follows 
that Y belongs to the great circle S^. It is clear that in this case m = 2. (For 
instance, m > 2 contradicts Theorem 2 for n = 2.) 

To prove our main results in this section for n = 3, 4 we need the following 
fact. (For n = 3, when A is an arc, a proof of this claim is trivial.) 



For m 




(5.1) 



A2 = 2/12/2 is an arc of length ijj. 



(5.2) 
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Lemma 2. Consider in S"~^ an arc lu and a regular simplex A, both are with 
edge lengths Tp, tp ^ it/2. Suppose the intersection of oj and A is not empty. 
Then at least one of the distances between vertices of u) and A is less than tp. 

Proof. We have u = U1M2, A = V1V2 ■ . - Vk, dist(Mi, U2) = dist{vi,Vj) = ^/j. 

Assume the eonverse. Then dist(?i,;, Wj) > ip for aU i, j. By U denote the 
union of the spherical caps of centers vi, i = 1, . . . , k, and radius ijj. Let B be 
the boundary of U. Note that ui and U2 don't He inside U. If {w/^, wJ,} = Ljf]B, 
then Tp = dist(ui, 7/2) > dist(?i'j^, (/J,), and P| A 7^ 0, where lu' = u'^ui^. 

We have the following optimization problem: to find an arc W1W2 of minimal 
length subject to the constraints wi,W2 S B, and W1W2 H A ^ 0? It is not hard 
to prove that dist(u>i,W2) attains its minimum when wi and W2 arc at the 
distance of ip from all Vi, i.e. WiVi . . .v^ and W2V1 . . .V/. are regular simplices 
with the common facet A. Using this, it can be shown by direct calculation 
that 

2kz'^ — (k — l)z — 1 

cosa = — — — , a = mmdist{wi,W2), z = costp (5.3) 

1 ~|~ [k I- j-^ 

We have a < tp. From (5.3) follows that cos a > z if and only if ^; > 1 or 
(fc + l)z + 1 < 0. It contradicts the assumption < z <1. □ 

5-C. Irreducible sets in S^. Now we consider irreducible sets for n = 3. In 
this case dim A^ < 2. 

Theorem 4. Suppose Y is irreducible and dim(Am) = 2. 

Then 3 < m < 5, and Am is a spherical regular triangle, rhomb, or equilateral 

pentagon with edge lengths tp. 

Proof. From Corollary 2 follows that ni < 5. On the other hand, m > 2. Then 

m = 3, 4, 5. Theorem 2 implies that A„i is a convex polygon with vertices 
2/1, ... , y„i- From Lemma 1 it follows that eo G Am, and degt/i ^ 1. 

First let us prove that if degt/i > 2 for all i, then Am is an equilateral m-gon 
with edge lengths ip. Indeed, it is clear for to = 3. 

Lemma 2 implies that two diagonals of A„i of lengths tp do not intersect each 
other. That yields the proof for m = 4. When to, = 5, it remains to consider the 
case where A5 consists of two regular non overlapping triangles with a common 
vertex (Fig. 3). This case contradicts the convexity of A5. Indeed, since the 
angular sum in spherical triangle is strictly greater than 180° and a larger side of 
spherical triangle subtends opposite large angle, we have Zyiyiyj > 60°. Then 

180° > Zj/2?/i2/5 = ^2/22/12/3 + ^ysyiVA + ^2/42/12/5 > 180° 
- a contradiction. 

Now we prove that degj/^ > 2. Suppose degyi = 1, i.e. (f>i,2 = ip, cpi^i > tp 
for i = 3,..., TO. (Recall that (pij = dist{yi,yj).) If eo ^ 2/i2/2j then after 
sufficiently small turn of 2/1 round 2/2 to Cq (Fig. 4) the distance 61 decreases - 
a contradiction. (This turn will be considered in Lemma 3 with more details.) 
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Fig. 3 Fig. 4 

It remains to consider the case: eo G 2/i2/2- If <Pi,j = where i > 2 or 
j > 2, then eo ^ ViVj- Indeed, in the converse case, we have two intersecting 
diagonals of lengths ^. Therefore, degj/i > 2 for 2 < i < m. For to = 3,4 it 
implies the proof. For m = 5 there is the case where Q3 = y^yiy^ is a regular 
triangle of side length ijj. Note that yiy2 cannot intersect (otherwise we 
again have intersecting diagonals of lengths V), then 2/12/2 is a side of A5. In 
this CclSCj cLS above, after sufficiently small turn of round j/2 to eo the distance 
di, i = 3, 4, 5, decreases - a contradiction. □ 

5-D. Rotations and irreducible sets in n dimensions. Now we extend 

these results to n dimensions.^ 

Let us consider a rotation R{ip, fi) on S"~^ about an (n — 3) - dimensional 
great sphere Vt in S"~^. Without loss of generality, we may assume that 

0. = {u= (wi, . . . , w„) e R" : ?ii = U2 = 0, Ml + . . . + = 1}. 

Denote by R{ip, fi) the rotation in the plane {wj = 0, i = 3, . . . , n} through an 
angle ip about the origin Q, : 

u'l = ui cos (fi — U2 sin (p, u'2 = ui sin tp + U2 cos tp, u'^ = Ui, i = 3, . . . ,n. 

Let 

H+ = {ue S"-^ : U2 > 0}, H_ = {u€ S""^ : U2 < 0}, 
Q = {uG S"-^ : U2 = 0, Ml > 0}, Q = {uG S"-^ : M2 = 0, Ml > 0}. 

^In the first version of this paper for m > n has been claimed that any vertex of T^{Y) 
has degree at least n — 1. However, E. Bannai, M. Tagami, and referees of this paper found 
some gaps in our exposition. Most of them are related to "degenerated" configurations. In 
this paper we need only the case n = 4,m < 6. For this case Bannai and Tagami verified 
each step of our proof, considered all "degenerated" configurations, and finally gave clean and 
detailed proof (see E. Bannai and M. Tagami: On optimal sets in Musin's paper "The kissing 
number in four dimensions" in the Proceedings of the COE Workshop on Sphere Packings, 
November 1-5, 2004, in Pukuoka Japan). Now this claim for all n can be considered only as 
conjecture. In 5-D we prove the claim when {j/j} axe in "general position". I wish to thank 
Eiichi Bannai, Makoto Tagami, and anonymous referees for helpful and useful comments. 
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Note that and arc closed hemispheres of S" ^, Q = QIJO, and Q is 
a hemisphere of the unit sphere 0,2 = {u ^ S"^^ : U2 = 0} boimdcd by O. 

Lemma 3. Consider two points y and eo in S"~^. Suppose y & Q and cq ^ Q. 
If Co G H^, then any rotation R{^p, fi) of y with sufficiently small positive {p 
decreases the distance between y and cq. 

If eo e H-, then any rotation R{(p,i}) of y with sufficiently small negative ip 
decreases the distance between y and cq- 

Proof. Let y be rotated into the point y{(p)- If the coordinate expressions of y 
and eo are 

y = (mi,0,M3, . . . Ml > 0; eo = (wi, W2, ■ • ■ , fn), then 

r{ip) := y{(p) ■ eo ~ uivi cos ip + uiV2 simp + uaus + . . . + UnVn. 
Therefore, r' {ip) = —uiVi sinip + uiV2C0S(p, i.e. r'(0)=uiV2- Then 

r'(0) > iff 1)2 > 0, i.e. eo GH+; 

r'(0) < iff V2<0 i.e. eo GH- ■ 

That proves the lemma for V2 ^ 0. In the case V2 ~ 0, by assumption (eo ^ Q) 
we have Ui < 0. In this case r'(0) = 0, and r"(0) = —UiVi > 0, i.e. = is a 
minimum point. This completes the proof. □ 

Proposition 2. Let Y be irreducible and m = |y| > n. Suppose there are no 
closed great hemispheres Q in S"~^ such that Q contains n — 1 points from Y 
and eo- Then any vertex ofT,^{Y) has degree at least n—1. 

Proof. Without loss of generality, we may assume that 

4>i,i=i^, i = 2, . . . ,degyi + 1; cpi^ > ip, z = degyi + 2, . . . , m. 

Suppose dcgyi < n — 1. Then > ■0 for i = n, . . . , m. Let us con- 
sider the great {n — 3) - dimensional sphere Cl in S"~^ that contains the points 
1/2, ■ . • , yn-i- Then Lemma 3 implies that a rotation R{>p, ft) of yi with suffi- 
ciently small ip decreases 6i. It contradicts the irreducibility of Y. □ 

Proposition 3. IfY is irreducible, \Y\ — n, dimA„ = n—\, then dcgyi = n—1 
for all i — 1, . . . ,n. In other words, A„ is a regular simplex of edge lengths ip. 

Proof. Clearly, A„ is a spherical simplex. Denote by Fj its facets, 
Fi := conv{?/i, . . . . . . ,2/„}. 

Let for cr C /„ := {1, . . . , n} 

Fa ■.= C]Fi. 
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We claim for j that: 



If eo ^ F{i,j}, then (f)ij = ip. (5.4) 

Conversely, from Lemma 3 follows that there exists a rotation R{(p,^ij) 
of yi (or yj if eo € Fj) decreases (respectively, 6j), where fly is the great 
(n — 3) — dimensional sphere contains F^^ jy . It contradicts the irreducibility 
assumption for Y. 

This yields, if there is no pair such that eo € then = V' for 

all i, j from 7„. 

Suppose eo G F„, where a has maximal size and |c7| > 1. Let a = In \ a. 
From (5.4) it follows that (pij = ijj if i € or j € a. It remains to prove that 
= tp for i,j e cr. 

Let A be the intersection of the spheres of centers yi, i G a, and radius tp. 
Then A is a sphere in S"~^ of dimension |(t| — 1. Note that =convex hull of 
{yi : i & a}, and for any fixed point x from Fa (in particular for x = eo) the 
distance dist(x,y) posses the same value (depending only on x) on the entire 
set y G A. TliCTi jji. i E a, lie in A at the same distance from cq. It is clciar that 

Y is irreducible if and only ii yi, i £ a, in A are vertices of a regular simplex of 
edge length tjj. 

Finally, we have that all edges of A„ are of lengths tp as required. □ 
Corollary 4. If n> 3, then A4 is a regular tetrahedron of edge lengths tp. 

Proof. Let us show that dimA4 = 3. In the converse case, dimA4 = 2, and 

from Theorem 4 follows that A4 is a rhomb. Suppose yij/3 is the minimal length 
diagonal of A4. Then (/)2,4 > tp (see Lemma 2). Let us consider a sufficiently 
small turn of the facet 2/12/22/3 round 2/12/3. If eo ^ yiy^, then this turn decreases 
either 9^ (if eo € 2/12/22/3) or ^2, a contradiction. In the case eo G 2/12/3 any turn 
of 2/2 round 2/12/3 decreases (p2,i and doesn't change ^2- Obviously, there is a 
turn such that (/i2,4 becomes is equal to tp. That contradicts the irreducibility of 

Y also. 

□ 

5-E. Irreducible sets in S^. 

Lemma 4. If Y cS^ is irreducible, \Y\ = 5, then degyt > 3 for all i. 

Proof. (1) Let us show that dim A5 = 3. In the converse case, dim A5 = 2, and 
from Theorem 4 follows that A5 is a convex equilateral pentagon. Suppose 2/12/3 
is the minimal length diagonal of A5. We have (p2,k > ■0 for fc > 3. Suppose 
eo ^ 2/12/3- If eo G 2/i2/22/3 then any sufficiently small turn of the facet 2/i2/32/42/5 
round 2/12/3 decreases 64 and ^5, otherwise it decreases 62, a contradiction. In 
the case Cq G yiys any turn of y2 round yiys decreases 02, fe for k = 4,5, and 
doesn't change 9i. It can be shown in the elementary way that there is a turn 
such that 02,4 or 02,5 becomes is equal to tp, a contradiction. 
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In three dimensions there exist only two combinatorial types of convex poly- 
topes with 5 vertices: (A) and (B) (see Fig. 5). In the case (A) the arc 2/3J/5 
lies inside A5, and for (B): 1/22/32/42/5 is a facet of A5. 




(A) (B) 

Fig. 5 

(2) By Sij we denote the arc yiyj, and by Sijk denote the triangle yiUjUk- Let Sijk 
be the intersection of the great 2— hemisphere Qijk and A5, where Qijk contains 
yi,yj,yk said boimdcd by the great circle passes through yi,yj. Proposition 2 
yields: if there are no k such that eo G s^fe, then degj/i > 3 for all i. 

It remains to consider all cases eo S Sijk ■ Note that for (A) Sijk ^ Sijk only 
for three cases: i = 1,2, 4; where j = 3, k = 5, or j = 5, k = 3 {siS5 = Sjss). 

(3) Lemma 1 yields that degyk > 0. Now we consider the cases degj/fe = 1, 2. 

If degyk = 1, (f>k,e = t/j, then eo e Ske- 

Indeed, otherwise there exists the great circle in such that fl contains 
yi, and the great sphere passes through fl and yk doesn't pass through eo. Then 
Lemma 3 implies that a rotation R{(f, O) of yk with sufficiently small (p decreases 
Ok - 8- contradiction. 

Since 9o < ip, cq can not be a vertex of A5. Therefore, eo lies inside Skt- 
Prom this follows if Sij for any j doesn't intersect Ske, then degj/j > 2. 

Arguing as above it is easy to prove that 

If degyk = 2, ^k,i = (t>k,j = tp, then eo e Sijk- 

(4) Now we prove that degyk > 2 for all k. Conversely, degyk = 1, eo € Ske- 

a) . First we consider the case when Ske is an "external" edge of A5. For the 
type (A) that means Ske differs from S35, and for (B) it is not S35 or S24- Since 
A5 is convex, there exists the great 2— sphere ^2 passes through yk, yi such that 
3 other points yi,yj,yq lie inside the hemisphere if+ bounded by fl2- Let fl be 
the great circle in O2 that contains ye and is orthogonal to the arc Ske- Then 
(Lemma 3) there exists a small turn of yi,yj,yq round O that simultaneously 
decreases 9i,9j,9q - a contradiction. 

b) . For the type (A) when degys = 1, 03,5 = ijj, eo G S35; we claim that S124 
is a regular triangle with side length tjj. Indeed, from a) follows that degy^ > 2 
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for i = 1,2,4. Moreover, if deg^j = 2, then cq = S35n'5i24- Therefore, in any 
case, (j)i^2 = 4>i,A = <t>2,A = We have the arc S35 and the regular triangle S124, 
both are with edge lengths V- Then from Lemma 2 follows that some < tp 

- a contradiction. 

c). Now for the type (B) consider the case: degj/3 = 1, ^3,5 = tp, cq G S35. 
Then for t/2 we have: degy2 = 1 only if (^2,4 = i>, eo = 524 0*35; deg2/2 = 2 
only if (f)2A = 4'2.b = ijj', and (f>2A = 4'i.2 = 4>2.b = ijj if dcgy2 = 3. Thns, in any 
case, (j)2,4 = ip- We have two intersecting diagonals 324,335 of lengths ^. Then 
Lemma 2 contradicts the assumption that F is a '^-code. This contradiction 
concludes the proof that degy/j > 2 for all k. 

(5) Finally let us prove that degyfe > 3 for all k. Assume the converse. Then 
degyk =2, eo G Sijk, where (j)k,i = <Pkj = fp- 

Case facet: Let Sijk be a facet of A5, and eo ^ Sij. By the same argument 
as in (4a), where CI2 be the great sphere contains Sijk, and CI be the great circle 
passes through yi,yj, we can prove that there exists a shift decreases 9e,9g for 
two other points ye, yg from Y, a contradiction. 

If eo G Sij, then any turn of sgq round Q doesn't change Og and 9q. However, 
if this turn is in a positive direction, then it decreases (pk,e and (pk,q- Clearly, 
there exists a turn when (p^^e or (pk,q is equal to V - a contradiction. 

It remains to consider all cases where Sijk is not a facet. Namely, there are 
the following cases: S124, S135 (type (A) and type (B)), S234 (type (B)). 

Case S124: We have degyi = 2, (pi^2 = (pi,4 = tp, eo G Si24- Consider a 
small turn of 7/3 roimd S24 towards yi. If eo ^ S24, then this turn decreases ^3. 
Therefore, the irreducibility yields (p^^^ = tp. In the case eo G S24, ^3 = ^3, but 
01,3 decreases. It again implies ^3^5 = tp. Since S35 cannot intersects a regular 
triangle S124 [see Lemma 2, (4b)], (p2A > ^- Then degy2 = degy^ = 3. (Since 
Co G S124, degy2 = 2 only if (p2,4 = tp.) Thus we have three isosceles triangles 
S243, S241, S245- Using this and ^3,5 = tp, we obviously have < for i = 3, 5, 

- a contradiction. 

Case Si35(type (B)) is equivalent to the Case S124. 

Ccise Si35(type (A)): This case has two subcases: S351, S153. 

In the subcase ,5i35 wc have dcgyi = 2, (^1.3 = 01,5 = tp, eo G S135. 
If eo ^ S135, then any turn of yi round .S35 decreases 9i (Lemma 3). Then 
eo G S135. Clearly, any small turn of 2/2 round S35 increases 02,4- On the other 
hand, this turn decreases 62 (if eo ^ S35) and (pi,2- Arguing as above, we get 
a contradiction. The subcase S315, where 03, 5 = tp, can be proven by the same 
arguments as Case si24- 

Case S234(type (B)): This case has two subcases: S243, S234- 

It is not hard to see that S243 follows from the Case facet, and S234 can be 
proven in the same way as the subcase §135. This concludes the proof. □ 

Lemma 4 yields that the degree of any vertex of r^(y) is not less than 
3. This implies that at least one vertex of T-^,{Y) has degree 4. Indeed, if all 
vertices of r^(y) are of degree 3, then the sum of the degrees equals 15, i.e. is 
not an even number. There exists only one type of {¥) with these conditions 
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(Fig. 6). The lengths of all edges of A5 except 3/22/4, 2/32/5 a^re equal to ■)/'• For 
fixed (/)2,4 — a, A5 is uniquely defined up to isometry. Therefore, we have the 
1-parametric family Ps^a) on S'^. If (^3,5 > </>2,4, then z > cosa > 2z — 1. 




Fig. 6: P5{a) 

Thus Theorem 4, Corollary 4 and Lemma 4 for n = 4 yield: 

Theorem 5. Let F C S'^ fee an irreducible set, \Y\ — ni < 5. Then Am for 
2 < m < A is a regular simplex of edge lengths ^, and A5 is isometric to -P5(a) 
for some a £ [V', arccos (2z — 1)]. 

5-F. Optimization problem. We see that if Y is optimal, then for some 
cases it can be defined up to isometry. For fixed yi e S"^^, « = 1, . . . ,m; the 
function H depends only on a position y ~ — j/o = eo G S"^^. Let 

H,n{y) /(I) + f{-y • 2/1) + • ■ • + f{-y ■ 2/m), 

i.e. Hm{y) ^ H{-y;Y). 

Thus for km we have the following (n— l)-dimcnsional optimization problem: 

km = max{i7™(?/)} 

V 

subjects to the constraint 

y e T{Y, 9o) := {y G A„ C S"-i : y ■ y, > to, i ^ 1, . . . ,m}. 
We present an efficient numerical method for this problem in the next section. 

6 On calculations of hm 

In this technical section we explain how to find an upper bound on km for 
n = 4, m < 6. Note that Theorem 5 gets for computation of hm a low- 
dimensional optimization problem (see 5-F) . Our first approach for this problem 
was to apply numerical methods |2S|. However, that is a nonconvex constrained 
optimization problem. In this case, the Nelder-Mead simplex method and other 
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local improvements methods cannot guarantee finding a global optimum. It's 
possible (using estimations of derivatives) to organize computational process 
in such way that it gives a global optimum. However, such solutions are very 
hard to verify and some mathematicians don't accept that kind of proofs. Fortu- 
nately, using geometric approach, estimations of hm can be reduced to relatively 
simple computations. 

Throughout this section we use the function f{9) defined for / e ^*{z) by 

^,„. _ J /(— 008 6*) < 6^ < ^0 = arccosfo (see Definition 4) 

Since / e ^*{z), f{6) is a monotone decreasing function in 9 on [0,^o]- 

6-A. The case m=2. Suppose m = 2 and Y is optimal for / G <P*(z). Then 
A2 = 2/12/2 is an arc of length tp, eo € A2, and 9i + 62 = V'l where 9i < 9o (see 
Lemma 1 and (5.2)). The efficient function F(6»i,6»2) = /(I) + f{9i) + /(6I2) is 
a symmetric function in 6'i , ^2 • 

We can assume that 9i < 92, then 9i e [ip — 9o,il)/2\. Since 62(^1) := tl) — 9i 
is a monotone decreasing function, /(O2(0i)) is a monotone increasing function 
in 9i. Thus for any 9i G [u, v] C [V' — V'/2] we have 

i=^(^i, ^2) < $2([«, v]):=f{l) + f{u) + /(^ - v). 

Let ui = ip — 9o, U2, ■ ■ ■ , UN, un+1 = ^/'/2 be points in [ip — 9o, ^/2] 
such that u,+i = Ui + s, where s = {9o —tjj/2)/N. If 9i G [u,,u,+i], then 
/i2 = H{yo;Y) = F{9r,92) < $2([«i,«i+i])- Thus 

h2 < X2{N,ip,9o) := max {$2(si)}, where Si := [ui,Ui+i]. 

l<i<N 

Clearly, X2{N,t{j,9o) tends to /12 as iV ^ 00 (e ^ 0). 

That implies a very simple method for calculation of /i2- Now we extend this 
approach to higher m. 

6-B. The function Qk- Suppose we know (up to isometry) optimal Y = 
{yi, ■ ■ ■ ,ym} C S"~^. Let us assume that dimA^ = n — 1, and V := convex 
hull of {yi . . . yn-i} is a facet of A^. Then rank{t/i, . . . , yn-i} = n — 1, and Y 
belongs to the hemisphere where contains Y and bounded by the great 
sphere S passes through V. 

Let us show that any y = y_^_ £ is uniquely determined by the set 
of distances 9i = dist{y, yi), i = 1, ... ,n — 1. Indeed, there are at most two 
solutions: y-^. G and y_ G H- of the quadratic equation 

y-y=l with y-yi = cos9i, i = 1, . . . ,n - 1. (6.1) 

Note that ?/+ = y_ if and only ii y E S. 

This implies that 9^, k > n is determined by 9i, i = 1, . . . ,n — 1; 

9k = &k{9i, ■ ■ ■ , 9n-l)- 
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It is not hard to solve (6.1) and, therefore, to give an cxphcit expression for Ok- 
For instance, let A„ be a regular simplex of edge lengths n/S. (We need this 
case for n = 3, 4.) Then ^ 

cos 6*3 = cose:,((9i,6l2) = 
= ^ ^cos 6i + cos 02 + y/Q — 8[cos^i cos ^2 + (cos ^2 — cos^i)^]^ ; 

cos (^4 — cos04(6'i, 6*2, ^3) — ^ ^cos6'i + COS 02 + cos 03+ 

10^ 1 + COS 9l COS 62 + COS 61 cos ^3 + cos 62 COS ^3 ~ ^ (cos^ 61 + cos^ 62 + cos^ ^3)^ . 

6-C. Extremal points of 9*; on D. Let a = (ai, . . . , a„_i), where < 
di < Oo < ip. (Recall that = dist{yi,yj); costjj = z; cos^o = to-) Now we 
consider a domain -D(a) in where 

Z)(a) = {y € H+ : dist(2/, yi) < ai, l<i<n - 1}. 

In other words, -D(a) is the intersection of the closed caps Cap(yi, ai) in H+: 

n-l 

£)(a)= nCap(yi,ai)f|i?+. 

Suppose dim£)(a) = n—1. Then -D(a) has "vertices", "edges", and "/c-faces" 
for < n — 1. Indeed, let 

a C I :={!,... ,n-l}, < |ct| < n - 1; 

Fa ■■= {y G D{a) : dist{y,y^) = W i€ a}. 

It is easy to prove that dimF^r = n — 1 — \a\; belongs to the boundary 
B of -D(a); and if cr C a', then F^' C F^. Actually, D{a) is combinatorially 
equivalent to an (n — l)-dimensional simplex. 

Now we consider the minimum of Qk{di, ■ ■ ■ ,dn-i) on £'(a) for k > n. In 
other words, we are looking for a point Pfc(a) g -D(a) such that 

dist(i/fc,pfc(a)) = dist{yk,D{a)). 

Since (fn^k > "4- > ^o- all Vk li(3 outside D{a). Clearly, 6^ achieves its minimum 
at some point in B. Therefore, there is ct C / such that 

Pfe(a) e Fa (6.2) 

Suppose <T = I, then F^ is a vertex of D{a). Let us denote this point by 
(a) . Note that the function 6/c at the point (a) is equal to (a) • 

*I am very grateful to referees for these explicit formulas. 
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Let crfc(a) denote cr C / of the maximal size such that cr satisfies (6.2). Then 
for cr/c(a) = /, Pfe(a) = p*(a), and for |crfe(a)| < n — 1, p/c(a) belongs to the 
open part of (a) • 

Consider n = 3. There are two cases for Pfe(a) (see Fig. 7): p?,{a) = p*(a) = 

2}) and P4(a) is the intersection in of the great circle passes through 
yi, 2/4, and the circle S{yi,ai) of center yi and radius ai (^'{i} C S{yi,ai)). 
The same holds for all dimensions. 

Denote by S'cr(fe) the great | a |— dimensional sphere passes through j/j, i £ a, 
and yk- Let <S'(yi, aj) be the sphere of center y^ and radius a,; and for cr C / 

So- := Pi S{yi,ai). 

Denote by s{a, k) the intersection of So-{k) and in 

s(a,/c)=5<,(A;)f|^<,f|ir+ 




Fig. 7 Fig. 8 



Lemma 5. Suppose £'(a) ^0, < < /o?" all i, and k > n. Then 

(i) Pk{a) e s(a-fe(a),fc), 

(m) if s{a, k) ^0, |0-| < n — 1, iften s(o-, fc) consists of the one point Pfe(a). 

Proof. («) Let := 6fe(pfe(a)) = dist(yfe,pfc(a)). Since 9fc achieves its minimum 
at J3fc(a), the sphere S{yk, 01) touches the sphere /S'cr(a) at Pfe(a). If some sphere 
touches the intersections of spheres, then the touching point belongs to the great 
sphere passes through the centers of these spheres. Thus Pfe(a) G So(a){k). 

[ii) Note that s{(T,k) belongs to the intersection in of the spheres S{yi,ai), 
i G a, and Scy{k). Any intersection of spheres is also a sphere. Since 

dimS'cr(fc) + dim^CT = n — 1, 

this intersection is empty, or is a 0— dimensional sphere (i.e. 2-points set). In the 
last case, one point lies in and another one in Therefore, s{a, k) = 0, 
or s{a, k) = {p}. Denote by a' the maximal size a' a such that s(cr'. A:) = {p\. 
It is not hard to see that S{^k, dist(yfc,p)) touches S^i at p. Thus p = pfe(a). □ 
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Lemma 5 implies a simple method for caknilations of the minimum of Ofe 
on -D(a). For this we can consider s{a,k), a <Z I, and if s{a,k) ^ 0, then 
s{a,k) = {pfe(a)}, so then 6^ attains its minimum at this point. In the case 
when An is a simplex we can find the minimum by very simple method. 

Corollary 5. Suppose \Y\ =n, < a, < for all i, and D{a) lies inside A„. 
Then 

On > 0n(ai, ■ • ■ , ffln-i) for all y G D(a). 

Proof. Clearly, A„ is a simplex. Since D{a) lies inside A„, for |a| < n — 1 the 
intersection of and Scr{k) is empty. Thus Pn(a) = p*(a). □ 

6-D. Upper bounds on Hm- Suppose dim = n — 1, and yi . . . j/„_i is 
a facet of A^- Then (see 5-F for the definition of Hm and T{Y, 6q)) 

Hm{y) = F{&1, ■ ■ ■ , ^n-l, ©n, • • • , ©m) = -^771(^*1, • • • , &n-l), 

where 

F„,{e,, On-l) ■= /(I) + M) + . . . + fiOn-l) + fiOniOl, • • • , ^n-l)) + • • • 

+ /(e„(^i,...,e„_i)). 

Lemma 6. Suppose f G ^*{z), \Y\ = m, dim A^ = n— 1, yi . . .yn-i is a facet 
of A„, dist(yi, yj) > > 9o for i^ j, < bi < Oi < 60 for i = 1, . . . ,n - 1; 
and efe(a) < 6*0 for all k > n. If £)(a) 0, then 

Hm{y)<^Y{h,a) for any y e E{h,a) := D{a)\U{h), 

where 

$K(b, a) := /(I) + fih) + ... + f{hn-i) + /(e„(p„(a))) + . . . + /(e„(p„(a))), 

n-l 

U{h) := y Cap(t/i,6i)- 

1=1 

Proof We have for 1 < i < n - 1 and y G E{h,a) : 9i > h (Fig. 8). By 
the monotonicity assumption this implies f{Oi) < fipi). On the other hand, 
y e D{a). Then Lemma 5 yields /(6ife) < f{@k{j>k{aS))) for k>n. □ 

Prom Corollary 5 and Lemma 6 follow 

Corollary 6. Let \Y\ = n. Suppose f, a, b, and Y satisfy the assumptions of 
Lemma 6 and Corollary 5. Then for any y e E{h, a) : 

H„,{y) < /(I) + /(61) + . . . + /(6„_i) + /(e„(a)). 
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Let K{n,9o) [0, 0o]"~'^i i-c K{n,do) is an (n — 1)— dimensional eubc 
of side length ^o- Consider for K{n, 9o) the cubic grid L{N) of sidelength e, 
where e = 9d/N for given positive integer N. Then the grid (tessellation) L{N) 
consists of A^"~^ cells, any cell c G L{N) is an {n — 1)— dimensional cube of 
sidelength e, and for any point (^^i, . . . , 9n~i) in c we have 

bi{c) <9i<ai{c), ai(c) = fei(c) +£, i = l,...,n-l. 

Let L{N) be the subset of cells c in L{N) such that D(a(c)) ^ 0. There exists 
c e -^(-^) such that if^ attains its maximum on T(Y,9o) at some point in 
E{h{c),a.{c)). Therefore, Lemma 6 yields 

Lemma 7. Suppose f and Y satisfy the assumptions of Lemma 6, N is a 
positive integer, and y G Am is such that dist(t/, yi) < 9o for all i. Then 

Hm{y)< inax {$r(b(c),a(c))} 

cez,(JV) 

6-E. Upper bounds on km- Suppose is a regular simplex of edge length ip. 
Then the efficient function F is a symmetric function in the variables 9\,. . . ,9m- 
Consider this problem only on the domain 

A := {y e A™ : V - < < ^2 < ■ • ■ < < ^o}- 

Let Lt^{N) be the subset of cells c in L{N) such that cf] A 7^ 0. If c e Lk{N) 
is such that i?(b(c), a(c))) lies inside A^,* then we have an explicit expression 
for $m(c) := $y (b(c), a(c)) (see Corollary 6). For n — A, Theorem 5 implies 
that Am is a regular simplex, where m = 2, 3, 4. Thus from Lemma 7 follows 

hm< >^m{N,il),6o) := max {$„(c)}. 

cGLa{N) 

Now we consider the case n = 4, m — b. Theorem 5 yields: A5 is isometric 
to Pnicx) for some a E '■— arccos (2z — 1)] (see Fig. 6). Let the ver- 

tices yi,y2,y3 of P5{a) be fixed. Then the vertices 1/4(0), 2/5(0) are uniquely 
determined by a. 

Note that for any y e D{9o,9o,9o) the distance ^4(0) := dist(t/, 2/4(0)) in- 
creases, and ^5(0) decreases whenever a increases. Let ai = V", 02, • ■ • , ctN, cun+i = 
tp' be points in [^p, tp'] such that a^+i = + e, where e = {tp' — ip)/N. Then 

6'4(ai) < 6'4(Q:i+i), 6*5(01) > 6'5(ai+i), 

so then 

/>4(o.)) > />4(o.+ i)), />5(o0) < />5(Oi+l)). 

Combining this with Lemma 7, we get 
/i5<A5(iV,V,^o) :=/(!) + 

max 2,3(c) + max {i?4 5(0, i)}}, 

ceL{N) l<i<N 
* Clearly, it holds for all c S La (AT) if AT is sufficiently large. 
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i?l,2,3(c) = /(61(C)) + /(62(C)) + /(63(C)), 

R4,5{c,i) = /(94(p4(a(c),ai))) + 7(65(^5(3(0), aj+i))), 

where pfe(a, a) = pfc(a) with yk = 2/fc(a). 

Clearly, Xm{^ + 1, V'l ^0) < ^mi^, V'l ^o)- It's not hard to show that 

/im = Ato(V',^o) := lim A„(iV,V,^o)- 

JV— >oo 

Finally let us consider the case: n = 4, m = 6. In this case, we give an upper 
bound on he by separate argument. 

Lemma 8. Let n = 4, / e ??*(2;), \/^ > > 2;, ^0 ^ [arccos y^, ^0] • Then 
he < max{/(e(,) + A5(V,^o), /(-\/i) + A5(V, ^0) }• 

Proof. Let F = {yi, . . . , yg} C C(eo, ^o) S^, where F is an optimal z-code. 
We may assume that ^1 < ^2 < • • • < ^6- Then from Corollary 3{i) follows that 

Oo> Oe> 0^> arccos-\/z. 

Let us consider two cases: (a) 9o > Oe > 0q, (b) 9q> Oe> arccos-^i. 

(a) We have he = H{yo;yi, . . . ,ye) = H{yo;yi, . . . ,^5) + /(^e), 

Hiyo; yi,..., 2/5) < /15 = A5(V, ^0), f{Oe) < f{e'o). 
Then /i6</>o)+A5(V',^o)- 

(b) In this case all 9i < 9q, i.e. Y c C(eo,^o)- Since 

H{yo;yu...,y5) < Mii'X), M) < /(-\/i), 

it follows that he < f{-y/z) + X^i'ip, 9'o). □ 

We have proved the following theorem. 

Theorem 6. Suppose n = 4, / e $*{z), ^/z > to > z > 0, and N is a positive 
integer. Then 

(i) ho = f{l), /ii = /(l) + /(-l); 

{a) hm = Xm{tp, 9o) < Xm{N, tp, 9q) for 2 < m < 5; 

(m) he < max{/>^) + A5(V',^o), /(-\/i) + A5(V',^o)} V e [arccos ^i, ^o] • 

6-F. Proof of Lemma B. First we show that fi e <?*(l/2) (see Fig. 9). 
Indeed, the polynomial /4 has two roots on [—1, 1]: ti = —to, to « 0.60794, t2 = 
1/2; /4(i) < for t £ [—to, 1/2], and /4 is a monotone decreasing function on 
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the interval [—1, —to] - The last property holds because there are no zeros of the 
derivative fi{t) on [-1,-to]. Thus, fi G $*{l/2). 

We have to > 0.6058. Then Corollary 3{ii) gives /U < 6. For calculations of 
hm let US apply Theorem 6 with tp = arccos2; = 60°, = arccosfo « 52.5588° 
We get 

ho = /(I) = 18.774, hi = /(I) + /(-I) = 24.48. 
H2 achieves its maximum at 61 = 30°. Then 

h2 = /(I) + 2/(- cos 30°) « 24.8644. 

For m = 3 we have 

/13 = A3(60°,6'o) « 24.8345 

at 63 = 00, 61=62!=^ 30.0715°. 

The polynomial H^i attains its maximum 

hi w 24.818 

at the point with 61 = 62^ 30.2310°, 6*3 = 6*4 « 51.6765°, and 

/15 « 24.6856 

at a = 60°, 61 « 42.1569°, 6^2 = 6*4 = 32.3025°, 6I3 = 6*5 = e^. 

Let O'q = 50°. We have /(50°) w 0.0906, arccos = 45°, /(45°) « 0.4533, 

A5(60°, Oq) = /15 « 24.6856, A5(60°, 50°) « 23.9181, 

/i6 < max{/(50°) + /15, /(45°) + A5(60°, 50°) } w 24.7762 < /12. 
Thus hmax = h2< 25. Since (4.2), we have S{X) < 25M. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



Fig. 9. The graph of the function f^it) 
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7 Concluding remarks 



This extension of the Delsarte method can be apphed to other dimensions and 
spherical ^-codes. The most interesting apphcation is a new proof for the 
Newton- Gregory problem, fc(3) < 13. In dimension three computations of hm 
are technically much more easier than for n = 4 (see |26p. 
Let 

^, , 2431 q 1287 7 18333 . 343 , 83 . 213 , t 1 

fit) = + + r H . 

■'^ ' 80 20 400 40 10 100 10 200 

Then / e to w 0.5907, /i(3,l/2,/) = 4, and h,nax = hi = 12.88. The 

(3) 

expansion of / in terms of Legendre polynomials Pk = G^. is 

/ = Po + 1.6Pi + 3.48P2 + 1.65P3 + 1.96P4 + O.IP5 + O.32P9. 

Since Cq = 1, Q > 0, we have fc(3) < h^ax — 12.88 < 13. 

Direct application of the method developed in this paper, presumably could 
lead to some improvements in the upper bounds on kissing numbers in dimen- 
sions 9, 10, 16, 17, 18 given in |21 Table 1.5]. ("Presumably" because the equality 
hmax — E is not proven yet.) 

In 9 and 10 dimensions Table 1.5 gives: 
306 < k{9) < 380, 500 < A:(10) < 595. 
Our method gives: 

n= 9 : deg/ = 11, P = /ii = 366.7822, to = 0.54; 
n = 10 : deg/ = 11, P = /ii = 570.5240, to = 0.586. 
For these dimensions there is a good chance to prove that 
fc(9) < 366, k{10) < 570. 

From the equality k{3) = 12 follows that <y93(13) < 60°. The method gives 
(/?3(13) < 59.4° (deg/ = 11). The lower bound on ipaiU) is 57.1367° 
Therefore, we have 57.1367° < ipailS) < 59.4°. 

Using our approach it can be proven that (^4(25) < 59.81°, (^4(24) < 60.5°. 
That improve the bounds: 

(^4(25) < 60.79°, (^4(24) < 61.65° 23 (cf. j4|); (^4(24) < 61.47° H; 

(y94(25) < 60.5°, (^4(24) < 61.41° 0. 
Now in these cases we have 

57.4988° < ipi{25) < 59.81°, 60° < (^4(24) < 60.5°. ^ 



However, for n = 5, 6, 7 direct use of this extension of the Delsarte method 
doesn't give better upper bounds on k{n) than Odlyzko-Sloane's bounds |77). 
It is an interesting problem to find better methods. 

^The long-standinding conjecture: The maximal kissing arrangment in four dimensions is 
unique up to isometry (in other words, is the "24-ceU"), and i/?4(24) = 60°. 
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Appendix. An algorithm 
for computation suitable polynomials f{t) 

In this Appendix is presented an algorithm for computation "optimal" ^ 
polynomials / such that f{t) is a monotone decreasing function on the interval 
[— 1, —to], and f{t) < for t £ [—to, z], to > z >Q. This algorithm based on 
our knowledge about optimal arrangement of points yi for given m. Coefficients 
Cfc can be found via discretization and linear programming; such method had 
been employed already by Odlyzko and Sloane for the same purpose. 

d , > 

Let us have a polynomial / represented in the form /(<) = 1 + ^kG^ {t). 

k=l 

We have the following constraints for /: (CI) Cfc > 0, 1 < k < d; 

(C2) /(a) > f{b) for -1 < a < 6 < -to; (C3) f{t) < for -to<t< z. 

We do not know cq where Hra attains its maximum, so for evaluation of hm 
let us use eg = j/c, where yc is the center of A^. All vertices y^ of A^, are at 
the distance of /9,„ from yc, where 



cos p„ 



= V(l + (m - 1)z)/t 



When TO = 2n — 2, Am presumably is a regular [n — l)-dimensional cross- 
polytopeJ In this case cospm — \fz- 

Let /„ = {1, . . . ,n} lj{2n - 2}, m e /„, 6m = -cospm, then 
Hm^Vc) = /(I) + If Fo is such that i/(yo; Y)<E^Fo + /(I), then 

(C4) ]{h^) <Folm, toG/„. Note that = Fq + 1 + ci + . . . + Cd = Fq + /(I) 
is a lower estimate of hmax- A polynomial / that satisfies (C1-C4) and gives 
the minimal E can be found by the following 

Algorithm. 

Input: n, z, to, d, N. 
Output: ci, . . . , Cd, Fq, E. 

First replace (C2) and (C3) by a finite set of inequalities at the points 
a,=-l + ej, 0<j<N, e=il + z)/N: 

Second use linear programming to find Ci, . . . , so as to minimize 

d 

E — 1 = Fo + J2 '^k subject to the constraints 
fe=i 



Ck>0, l<k<d- ^CfeGi"^(aj)>^CfeGi"^(aj+i), a, £ [-1, -to]; 



k=l k=l 
d 



H-EcfeGi"^(a,) <0, a,e[-to,z]; 1 + ^ CfeG^") (5„) < Fq/to, to g /„. 

k=l k=l 

Let us note again that E < h^ax, and E = h^nax only if h^ax = Hoodie) 
for some toq G In- 

^Open problem: is it true that for given io, d this algorithm defines / with minimal hmax'^ 
^It is also an open problem. 
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